
Finding Saddlepoints of Two-Person, Zero Sum Games 
Author(s): Donna Crystal Llewellyn, Craig Tovey, Michael Trick 

Source: The American Mathematical Monthly , Yol. 95, No. 10 (Dec., 1988), pp. 912-918 
Published by: Mathematical Association of America 
Stable URL: http://www.jstor.org/stable/2322384 

Accessed: 25/03/2009 14:31 


Your use of the JSTOR archive indicates your acceptance of JSTOR's Terms and Conditions of Use, available at 
http://www.jstor.org/page/info/about/policies/terms.jsp. JSTOR's Terms and Conditions of Use provides, in part, that unless 
you have obtained prior permission, you may not download an entire issue of a journal or multiple copies of articles, and you 
may use content in the JSTOR archive only for your personal, non-commercial use. 

Please contact the publisher regarding any further use of this work. Publisher contact information may be obtained at 
http ://www.j stor.org/action/sho wPublisher?publisherCode=maa. 

Each copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed 
page of such transmission. 

JSTOR is a not-for-profit organization founded in 1995 to build trusted digital archives for scholarship. We work with the 
scholarly community to preserve their work and the materials they rely upon, and to build a common research platform that 
promotes the discovery and use of these resources. For more information about JSTOR, please contact support@jstor.org. 


Mathematical Association of America is collaborating with JSTOR to digitize, preserve and extend access to 
The American Mathematical Monthly. 


STOR 


http ://www.j stor.org 





Finding Saddlepoints of Two-Person, Zero Sum Games 


Donna Crystal Llewellyn 1 and Craig Tovey 2 , 

School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta 


Michael Trick 3 , Carnegie-Mellon University 
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1. Introduction. A two-player, zero sum game, as defined by von Neumann and 
Morgenstem in their classical work [13], is completely specified by its payoff matrix. 
Let player I have m possible “moves” or “plays” (called pure strategies) and player 
II have n pure strategies. Suppose that player I pays player II a iJ if player I chooses 
his zth strategy and player II chooses her jth. Then the payoff matrix of the game is 
the m X n matrix A = [a t j\. 
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A saddlepoint (SP) of a matrix A is an entry a ij9 where a tj is the largest element 
in row / and the smallest in column j; that is 

a ij ^ a ik f° r ^ 
a iJ - < a kj for all k. 

If both inequalities are strict, then we call a ij a strict saddlepoint (strict SP). 

If there exists a strict saddlepoint of the payoff matrix, then it is the unique 
optimum solution to the game. A saddlepoint of the payoff matrix (if it exists) 
corresponds to a (not necessarily unique) optimum solution. Notice that if all of the 
entries of the payoff matrix are distinct, then a saddlepoint is a strict saddlepoint. 
(In general, finding the optimal solution vectors for the two players is equivalent to 
solving an m X n linear programming problem [3].) 

While these properties of saddlepoints are well known, it appears that the 
question of how to find these solutions efficiently has not been carefully studied. 
For instance, in [ 4 ] it is stated that “a [saddlepoint] is easy to find, even in a very 
large matrix.” Other standard texts, such as [ 1 ], [ 2 ], [ 5 ]— [ 10 ], [ 12 ] and [ 14 ]— [ 17 ] also 
mention that saddlepoints are easy to locate but fail to discuss finding an efficient 
algorithm. 

The purpose of this article is to show that a strict saddlepoint, when it exists, can 
be found very quickly. In fact, we give an algorithm that requires looking at only 
— w 0 58^ entries of the payoff matrix (where m > n). Hence as the number of 
strategies increases, the proportion of the payoff matrix examined approaches zero. 
This algorithm also determines when a matrix has no strict SP. Further, we will 
prove that finding a (nonstrict) saddlepoint requires examination of the entire 
payoff matrix. 

Given a matrix, we refer to the problem of finding a (strict) saddlepoint or 
determining that one does not exist as the (strict) SP problem. For ease of notation, 
we will denote log 2 x by lg x. 

We distinguish between two different types of computations. One kind is the 
ordinary arithmetic operation (i.e., comparison, addition, multiplication and divi- 
sion by integers). The other kind is the identification of the value of a matrix entry; 
this might involve a lengthy function evaluation. Therefore, we adopt the following 
model of computation. Assume that each matrix query contributes a units to the 
running time and each arithmetic operation requires ft time; where a and fi are 
some constants. 

II. Finding Strict Saddlepoints. We motivate our algorithm for finding a strict 
saddlepoint with a couple of observations. 

Observation 1. If two entries of A, a iJ and a kl are known, then at least one of them 
can be eliminated as a possible strict SP with a single matrix look-up and a constant 
number of arithmetic operations. 

Proof Clearly this follows if either i = k or j = /. Hence assume that neither of 
these holds. Without loss of generality let a tj < a kl . Look up a u . Clearly, if a ij is a 
strict SP then a tJ > a u . Further, if a kl is a strict SP then a H > a kl . Since a tj < a kh 
it is clear we can not have a tj > a a > a kl . Hence, knowing a a will eliminate either 
a or a kl or both from consideration. 
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Notice that in the proof of observation 1 it is assumed that a ij < a kl . Hence, the 
entry a kj cannot be a strict SP since a ij < a kl implies that it is impossible to have 
a tj > a kj > a kl . This leads to our second observation. 

Observation 2. If the diagonal terms of a square matrix, A, are known and are such 
that a u < a 22 < • • • < a nn , then no term below the diagonal can be a strict SP. 

Note that to apply observation 2 we will first have to rearrange the rows and 
columns of A so that the diagonal is ordered (which will not affect the existence or 
nonexistence of a strict SP). Our algorithm uses observations 1 and 2 together with 
proper manipulation of rearrangement pointers. The idea of the algorithm is 
illustrated in Figures 1 and 2. 



Fig. 1 


Fig. 2 


The main diagonal is looked up and then sorted to eliminate the subdiagonal 
portion of the matrix from consideration. The remaining region (i.e., the upper 
triangle) is divided into two triangles and one square which are studied recursively 
(see Figure 1 for the first step of this recursion and Figure 2 for the next two 
steps). These pieces are then put back together using observation 1. In the diagrams, 
the shaded regions are the parts of the matrix that have been ruled out by our 
observations and hence do not have to be searched. 

We now state the algorithm precisely. 


Procedure Square 

Input: A , n , r(l), . . . , r(n ), c(l), . . . , c(n ) 
Output: candidate 


This procedure takes as input a matrix A , and a subset of n of its rows and a subset 
of n of its columns. These rows and columns define_an n X n submatrix A of A. 
The procedure then returns candidate, a location in A which might be a strict SP of 
A. No other element of A is a strict SP of A. 

Step 0: If n = 1, return candidate = (r(l), c(l)) 

Step 1: Look up and sort A{r{i\ c(i)), i = in ascending order of size. 
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Suppose a is the permutation that gives this sorting, i.e., suppose that the set 
{ A ( r(o( i ))), c( a(i ))}"_ , is an increasing sequence. 

Set 

r(i) = r(a(i)) 

c(i) = c(o(i)) for / = 1, . . . , n. 

Step 2: Call Triangle (A, n , r(l), . . . , r(n ), c( 1), . . . , c(n); candidate) 

Step 3: If candidate = 0, return 0. Otherwise, take a -1 (candidate) and return it. 


Procedure Triangle 

Input: A , n , r( 1), . . . , r(n), c(l), . . . , c(n) 
Output: candidate 


This procedure takes as input a matrix A and a subset of n of its rows and a subset 
of n of its columns. These rows and columns define an upper triangle submatrix T 
of A. The procedure then returns a location in T which might be a strict SP of A. 
No other element of T is a strict SP of A. 

Step 0: If n = 1, return candidate = (r(l), c(l)) 

Step 1: Call Triangle (A, If J, r( 1), . . . , rflf J), c(l), . . . , c( If J); candidatel) 

Step 2: Call Square (A, [f ], r(l), . . . , r([f 1), c([f ]), . . . , c(«); candidate2) 

Step 3: Call Triangle (v4, [f J), r(l + ffl, , r(n ), c(l + ff ]), . . . , c(«); candidate3) 
Step 4: Using Observation 1, discard at least two of the candidates. If candidate/ 
remains (/ = 1, 2 or 3), let candidate equal candidate/. Otherwise, candidate = 0. 
Step 5: Return candidate. 

Program Strict Saddlepoint 

This program solves the strict saddlepoint problem on an n X n matrix A. 

Step 1: Call Square (n, 1, . . . , n, 1, . . . , n; candidate) 

Step 2: If candidate =£ 0, check to see if candidate is a strict SP; if so, then report 
the strict SP. Otherwise, then conclude there does not exist a strict SP. 

Proposition. Program Strict Saddlepoint solves the strict saddlepoint problem for 
a square matrix. 

Proof Observation 2 implies that candidatel, candidate2, and candidate3 are the 
only possible positions of a strict saddlepoint in the upper triangle matrix spanned 
by die row and column indices which, are passed to Procedure Triangle. Let 
candidate be the output from Procedure Square. Then by the above remark and 
Step 3 of Procedure Square, this position is the only possible position in the square 
matrix defined by the indices passed to this procedure by the central program. Thus 
all that is left for the Program Saddlepoint to do is to check if this position is indeed 
a strict saddlepoint; which it does in Step 2. 

Theorem. Program Strict Saddlepoint terminates in time 0(n {g 3 ). 

Proof Let F(n) be the time required to run Procedure Square 
(«, r(l), . . . , r(n ), c( 1), . . . , c(n); candidate) and let G(n) be the analogous time for 
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procedure Triangle (n, r(l), r(n), c(l), . .., c(n); candidate). Hence, the work 
required to run Program Saddlepoint is F(n) 4- O(fin); (due to the last check 
required to see if candidate is really a strict SP in step 2 of Program Saddlepoint). 
Notice that G(n ) = F(n) — c 1 /in log (/in) for some constant, c v because the only 
difference in the two algorithms is that the sorting step of Procedure Square does 
not occur in Procedure Triangle. By analyzing the steps of Program Saddlepoint, 

F( n ) = c ift n lg($ n ) + F ( 2) + 2 G (l) + + 

where the c 3 fin term is the cost of passing the row and column pointers and c 2 /i is 
the cost of making the final decision between candidates; that is, step 4 of the 
triangle procedure (using observation 1). 

So, 

F(n ) = c^n lg (/?«) + 3F^ -) + c 2 /J + c 3 /i« - 2c i^( 2 ) 18 (t~) 

= c x n lgn + 3 f| -) + c 2 0 + c 3 (0«) - 


= 3f(-J + c lP + (c, - Cj)(pn). 

Iterating this expression yields 

F(n) = 3 lg "F(l) + S £ ^ + (^) (c 3 - Cl )(/8«)J. 

But, F(l) = a, so after simplifying this geometric series we get 

F{n) = 3 Ig " | a + ^c 2 js) + 2^) 8 (c 3 - Cl )(pn) - {^c 2 ! 3 + 2 (c 3 - q)/?*) 

= 3 lg "|a + —c 2 /3 + -(c 3 - - | ~c 2 P + 2 (c 3 - c^nj. 

So, 

F(n) = 0( 3 lgn ). 

But, 3 lgn = n lg3 , and hence 


F(n) = 0(n lg3 ). 

Hence, by the remark at the beginning of the proof, the Program Strict Saddlepoint 
runs in time 0(n Xg3 ) as claimed. 


Note that this analysis shows that the program saddlepoint requires both order 
n i.5s look-ups and order n l 5S total running time. That is, even if a matrix 

look-up takes significantly more time than an arithmetic operation, that is, a » /}, 
the algorithm will remain 0(n 15S ). 

The reader familiar with fractals will note that the unshaded region in Figure 2, 
corresponding to the portion of the matrix that must be searched, is an illustration 
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of a fractal of dimension lg 3; and hence has area « lg3 . This provides a geometrically 
intuitive argument for the running time of Program Strict Saddlepoint [ 11 , pp. 
141-142]. 

Corollary. The strict SP problem on an m X n matrix , with m > n, can be 
solved in <9((™)/7 lg3 ) time. 

Proof. Divide the matrix into \f] possibly overlapping n X n square matrices, 
A v ..., A [m/n y Apply the theorem to each A i to get at most \f] candidate strict 
SP’s, from which the real strict SP (if it exists) can be found in work 0{~ + m + n). 
The total work is thus bounded by 


m 

n 


0(n lg3 ) + O(m) = O 



,lg3 


as claimed. 

III. Saddlepoints. Note that if we are seeking a SP then observation 1 holds only 
if a n > a 22 > • • • > a nn and similarly observation 2 will fail unless the inequality 
is strict. Thus it is not clear whether the search for a SP can be sped up in a similar 

way. In fact we show next that the SP problem cannot be solved in time better than 

2 

cn . 

Let C be the class of all algorithms that are made up entirely of table look-ups 
and arithmetic operations. 

Theorem. There does not exist an algorithm for the SP problem in the class C that 
can run faster than cn 2 . Further , there exists an algorithm in C that does run that fast. 

Proof We can put a check mark in each matrix cell which is the maximum in its 
row in time cn 2 , and then do the same for the cells which are the minimum in their 
columns in cn 2 time while checking for a cell with two check marks which would 
clearly be an SP. Hence, there surely is a quadratic time algorithm in C for this 
problem. 

We now show that no faster algorithm is possible. We do this by playing the 
adversary against an arbitrary algorithm in C. We argue that the algorithm must 
look up every entry in A. 

Our strategy: For the first n 2 - 1 queries of the values of a tj we answer 0, unless 
all the other entries in row i have been queried, in which case we answer 1. If we 
answer .5 to the last (n 2 ) query, then either that cell is the unique SP or there does 
not exist an SP. This follows since every 0 entry has either a 1 or a .5 in its row and 
every 1 has either a 0 or a .5 in its column. If instead, we answer — 1, then the cell is 
not' a SP, and in fact any zero entry in that last row will be a SP. Hence no 
algorithm in C can know the answer to the SP problem without n 2 queries. 

In conclusion, we have shown that the strict SP problem can be solved in 0(« L58 ) 
time, but that a similar improvement from 0(n 2 ) for the SP problem cannot be 
found. Note that our algorithm for the strict SP problem uses 0(« 158 ) matrix 
queries and has an 0(n l5S ) total running time. We conjecture that it may be 
possible to solve this problem with less than 0(n l 5S ) matrix queries but at the cost 
of a longer running time. 



918 


D. C. LLEWELLYN, C. TOVEY, AND M. TRICK 


REFERENCES 

1. Jean Pierre Aubin, Mathematical Methods of Game and Economic Theory, North- Holland Publish- 
ing Company, Amsterdam, 1979. 

2. David Harold Blackwell and Meyer A. Girshick, Theory of Games and Statistical Decisions, Wiley, 
New York, 1954. 

3. George B. Dantzig, Linear Programming and Extensions, Princeton University Press, Princeton, NJ, 
1963. 

4. Melvin Dresher, L. S. Shapley, and A. W. Tucker, editors, Advances in Game Theory, Princeton 
University Press, Princeton, NJ, 1964. 

5. David Gale, The Theory of Linear Economic Models, McGraw-Hill, New York, 1960. 

6. David Gale, The Theory of Matrix Games and Linear Economic Models, Department of Mathe- 
matics, Brown University, Providence, RI, 1957. 

7. Samuel Karlin, Mathematical Methods and Theory in Games, Programming, and Economics, 
Addison- Wesley Publishing Company, Reading, MA, 1959. 

8. H. W. Kuhn and A. W. Tucker, editors, Contributions to the Theory of Games, Volume I. Princeton 
University Press, Princeton, NJ, 1950. 

9. Richard I. Levin and Robert B. DesJardins, Theory of Games and Strategies, International 
Textbook Company, Scranton, PA, 1970. 

10. Duncan R. Luce and Howard Raiffa, Games and Decisions: Introduction and Critical Survey, John 
Wiley and Sons, Inc., New York, 1957. 

11. Benoit B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman and Company, New York, 
1977. 

12. Francis B. May, Introduction to Games of Strategy. Allyn and Bacon, Boston, MA, 1970. 

13. John von Neumann and Oskar Morgenstem, Theory of Games and Economic Behavior, Princeton 
University Press, Princeton, NJ, 1944. 

14. Guillermo Owen, Game Theory, second edition, Academic Press, New York, 1982. 

15. T. Parthasarathy and T. E. S. Raghavan, Some Topics in Two- Person Games, American Elsevier 
Publishing Company, New York, 1971. 

16. Anatol Rapoport, Two Person Game Theory; the Essential Ideas, University of Michigan Press, 
Ann Arbor, 1966. 

17. S. Vajda, The Theory of Games and Linear Programming, Methuen, London, 1967. 



